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Abstract 

We evaluate the squeezing of a probe beam with a transverse Gaus- 
sian structure interacting with an ensemble of two-level atoms in a cavity. 
We use the linear input-output formalism where the effect of atoms is 
described by susceptibility and noise functions, and show that the trans- 
verse structure is accounted for by averaging these atomic functions over 
the intensity profile. The results of the plane-wave and Gaussian-wave 
theories are compared. We find that, when a large squeezing is predicted, 
the prediction of the plane-wave model is not reliable outside the Kerr 
domain. We give an estimate of the squeezing degradation due to the 
Gaussian transverse structure. 

1 Introduction 

When an optical cavity containing a nonlinear medium approaches a bistability 
turning point [Q , the fluctuations of the incoming beam are squeezed Q . Nearly 
perfect squeezing is expected when the nonlinear medium is a lossless Kerr 
medium || Squeezing was experimentally demonstrated at the output of 
a bistable cavity containing an atomic beam (5[ [| or a cloud of cold atoms 
released from a magneto-optical trap (7). A lot of theoretical work has been 
devoted to the case where the nonlinear medium consists in two-level atoms 

Most theoretical treatments have simplified the description of the field- 
matter interaction by modelizing the laser beam as a plane-wave. It is how- 
ever clear that the description of the bistable behaviour of the cavity has to 
be modified to account for the transverse Gaussian profile of the laser beam 
[ fl7[ |l8 , 19, 20]. This is also true for the predictions concerning field fluctuations 



at the output of the cavity, as shown in references |2lj, |22|, |23j. These papers 
contain a derivation of squeezing spectra in optical bistability accounting for 
the transverse Gaussian structure of the laser beam. For particular sets of pa- 
rameters, the spectra obtained in Gaussian-wave calculations are compared by 
the authors with the predictions of the plane-wave model. In the good-cavity 
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limit where the atomic damping occurs on a much shorter time than the in- 
tracavity field damping and where the atomic variables can be adiabatically 
eliminated, Xiao et al |2lJ conclude that the plane-wave calculations may be 
trusted in either the weak- field or dispersive limits. In a complete treatment 
without adiabatic elimination , Hope et al relate squeezing degradation due 
to transverse structure to the ratio of atomic and cavity damping times. 

In the present paper, we quantitatively delineate the regions in the parameter 
space where the predictions of the plane-wave and Gaussian-wave calculations 
coincide, and therefore estimate the reliability of plane-wave model. We obtain 
criteria for reliability of plane- wave model which are quite different from those 
of references jil], ||[ . 

We will use the linear input-output formalism generalized to incorporate the 
treatment of atomic fluctuations []l6| . This formalism is based upon a linearisa- 
tion of the fluctuations around the working point of the bistable cavity, as the 
quantum stochastic methods [^4| , and it provides the same spectra as the latter 
when the same system is studied with the same simplifying assumptions p5[ . In 
contrast with the quantum stochastic methods, it provides an intuitive under- 
standing of the noise processing occuring in the bistable cavity. The squeezing 
process through the parametric transformation of vacuum fluctuations going 
into the intracavity mode is described in terms of atomic susceptibility func- 
tions, while the squeezing degradation through the addition of fluctuations by 
the atoms is characterized by atomic noise functions. 

We will show in the following that the effect of Gaussian profile upon noise 
processing appears quite naturally as a spatial averaging of these susceptibility 
and noise functions. This will allow us to obtain simple techniques for describing 
the influence of transverse structure. It is worth stressing that the property of 
averaging of the susceptibility and noise functions is not postulated a priori, 
but is rather a consequence of the simple assumptions made in the treatment of 
Gaussian profile: we assume that the higher-order transverse modes are either 
far from resonance with the cavity, or of such a high order that they are not 
efficiently coupled to the fundamental mode. It is thus possible to project any 
field radiated by the atoms upon this mode, and to obtain the averaging property 
as a consequence of this projection. 

Using this technique, we will first recover the well-known result 0, ^|, 
[l9L ||(| that the mean fields evaluated with the plane-wave and Gaussian-wave 
calculations coincide in the low-saturation limit where the population of the 
atomic excited state remains small. This implies that the bistability curves 
relating the input and intracavity intensities also coincide provided that these 
intensities are properly defined. We will then evaluate the atomic susceptibility 
and noise functions, and deduce that the conditions for coincidence of plane- 
wave and Gaussian-wave calculations are much more restrictive in this case: 
the limit of a small population in the excited state is not sufficient for ensuring 
the coincidence for atomic noise functions. In particular, we will conclude that 
the noise functions evaluated in the plane-wave model are not reliable when 



2 



large squeezing is predicted with parameters lying outside the Kerr domain. 
We will finally give a quantitative estimate for the degradation of squeezing due 
to transverse structure. 



2 Model and assumptions 

We consider a single ended ring cavity containing a nonlinear medium. A probe 
beam is entering the cavity through the coupling mirror, interacting with the 
medium and leaving the cavity through the same mirror. To study the influence 
of this device on the quantum fluctuations of the probe beam, we use a semi- 
classical input-output formalism, in which the field fluctuations are considered 
as classical stochastic variables, which are driven by the quantum fluctuations of 
the incoming fields at the laser frequency lul- The linear input-output relations 
give the quantum fluctuations of the output field as a function of those of the 
input field, in cases when they are small compared to the mean field-values, 
an assumption which is particularly well adapted for high photon numbers, i.e. 
laser beams in a resonator. The method is described in detail in Jig ], 

The intracavity field is defined as Heisenberg operator in the frame rotating 
with respect to the laser frequency lu^: 

e -i^E(r,t)+e iULt E(r,t)^ (1) 

In order to account for the spatial beam structure, the components E (r, t) 
and E{r,t)^ are written as products of a Gaussian mode function by a time- 
dependent mode amplitude (r stands for the three-dimensional position (x, y, z)): 



E(r,t) = x ^A(t)u(r) (2) 
V 2e c 

E(r,ty = ^A^uivf (3) 

This definition is such that (^A (i)^ A (t)^ is the number of photons going through 
a beam section in the unit time. 

The normalized Gaussian mode functions are written for a propagation along 
the x direction: 



U (r) = y^exp[-^(r)] (4) 

/(?') is the Gaussian intensity normalized to unity on the beam axis, S the 
effective beam section and f(r) the phase: 
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(5) 

(6) 
(7) 



w (x) denotes the position-dependent beam size, wq the beam waist, Ir the 
Rayleigh divergence length and A the laser wavelength: 



w 2 (x) = ™ 2 o[1 + Jt) (8) 

In = ( 9 ) 

Since we suppose the cavity to be perfectly mode matched to the laser beam, 
the same transverse structure holds for the input field and for the intracavity 
field before its interaction with the atoms. This structure is modified for the 
output field however, through the non linear interaction with the atoms on one 
hand, and through spontaneous emission of light in all possible directions by the 
atoms on the other hand. Both effects couple the fundamental Gaussian mode 
to all the other cavity modes. We will suppose here that the field radiated by the 
atoms into the other modes has an effect on intracavity fields much smaller than 
the field radiated into the fundamental Gaussian mode. This may be due either 
to the fact that the fundamental Gaussian mode is close to a cavity resonance 
while the other modes at the same frequency are far from any cavity resonance, 
or that they have such a high order that their coupling with the fundamental 
mode is negligeable. We will also consider that the detection setup is perfectly 
matched to the Gaussian mode. 

As a consequence, we can write the input and output fields E m (r,t) and 
E out (r,t) in the same manner as the intracavity field E(r,t): 



E ia (r,t) = ^L^{t)u{v) (10) 

£ out M = v1S A ° ut(<)u(r) 

Clearly the reflection-transmission equations for the fields on the coupling mirror 
will be the same for the field amplitudes as for the plane- wave fields in |l(J ■ In 
contrast, the effect of the atoms on these field amplitudes will be modified, due 
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to the non-linearity of the interaction and to the Gaussian variation of laser 
intensity. In order to evaluate the modification of the mode amplitudes due to 
the atomic medium, we will disregard the fields radiated by the atoms into the 
higher-order modes and obtain these modifications by projecting the radiated 
field onto the fundamental mode. This will be the key simplifying assumption 
in the forthcoming calculations concerning the mean fields as well as the field 
fluctuations. 

We will also use two simplifying hypothesis. First, we will consider that the 
transverse size w of the probe beam is much smaller than the size of the atomic 
medium, so that the dependence of the atomic density versus the transverse 
variables y and z may be disregarded. This density p (r) will therefore be 
replaced by the on- axis atomic density p m : 

p m (x)=p(x,y = 0,z = 0) (12) 

Note that the longitudinal dependence of the atomic density will have to be 
accounted for, since we consider that the length of the interaction zone is limited 
by the size of the atomic medium. Second, we will assume that the Rayleigh 
divergence length Ir is much larger than this size, so that the x-dependence of 
the laser beam size w can be ignored. However, we will not suppose that the 
atomic medium is centered on the cavity waist. 



3 Description of mean fields in a Gaussian-wave 
theory 

The atomic medium is considered to be a homogeneously broadened system of 
two-level atoms with a resonance frequency uj , which is placed inside an optical 
cavity and driven by a coherent mean field of frequency u>l close to a cavity 
mode u>c- Atom- laser detuning S and cavity-laser detuning <j) are respectively 
normalised by the atomic dipole decay rate 7 and the cavity field decay rate k: 



6 = (13) 

7 

4> = (14) 

It follows from the discussions of the previous section that the evolution of 
the intracavity Gaussian mode amplitude during one cavity round trip time is 
described by the same differential equation as in the plane-wave model: 



T C d t A(t) = -ktc (1 + i<f>) A(t) + V^cA m (t) + dA(t) (15) 

However, the modification dA of the intracavity Gaussian mode amplitude due 
to one passage through the atomic medium has not the same expression as 
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in the plane-wave model, as a consequence of the nonlinearity of atoms-field 
interaction. In the present section, we evaluate its mean value and then deduce 
the bistability curve. 

The effect of an optically thin atomic layer of length dx may be described 
by a local modification dE (r) of the field: 

E (r) -> E (r) + dE (r) (16) 

with f|: 

Q\2 

dE(r) =-— p(r)a(r)E(r)dx (17) 
An 

a (r) is the atomic polarizability derived from the solution of the optical Bloch 
equations fl6| and measured as a dimensionless number: 

on is the dimensionless linear polarizability of the atomic dipolc: 

« = ihs (19) 

s (r) is the dimensionless saturation parameter proportional to the local laser 
intensity: 



•W = ftF (20) 



p (r) is the dimensionless field parameter, that is precisely the Rabi frequency 
normalized to the decay rate of the dipole: 



„. . d EM /3A 2 A 

/3(r) - J ¥T = V^vf M(r) (21) 

where do is the matrix element of the atomic dipolc. These relations may also 
be written in terms of the saturation parameter s m and of the field parameter 
f3 m evaluated on the beam axis: 



»(r) = s m /(r) (22) 

- ^ rrt < 23 » 
«■ = < 24 » 

Note that we may forget the ^-dependence of the on-axis saturation parame- 
ter s m as a consequence of the assumption of a large Rayleigh divergence length. 
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The more general case of an arbitrary Rayleigh divergence length might as well 
be studied p7| . Using the transverse monomode approximation discussed above, 
we project dE onto the Gaussian mode, and write the following equation for 
the modification dA of the mode amplitude, summing up over the layers inside 
the atomic medium: 

dA 3A 2 

-A = ~A^ ains (25) 
n s is the effective number of atoms in the probe beam, taking into account the 
effect of the saturation: 

In order to evaluate n s , we expand it in powers of the saturation parameter: 



£(-2 Sm )V 1+fe ) (27) 

k=0 

J f{v) 1+k p{v)dv (28) 



nA 1 ' , that we will denote simply n, is the effective number of atoms in the probe 
beam, defined in the unsaturated regime: 

n = J f(v)p{v)dv (29) 

The next term corresponds to the lowest-order saturation effect, that is 
precisely the Kerr nonlincarity of the atomic medium. The other terms represent 
higher-order saturation effects. Due to the Gaussian shape of the function / (r), 
the function / (r) 1+ has the same expression as / (r) with a modified value of 
the beam size parameter, precisely with w 2 replaced by w 2 / (1 + k); n^ 1+fe ' has 
therefore the same expression as n with w 2 replaced by w 2 / (1 + k). 

We now use the assumption that the transverse size w of the probe beam 
is much smaller than the size of the atomic medium. The unsaturated effective 
number n has thus a simple expression in terms of the on-axis atomic density: 



irw 2 



p m (x) dx (30) 



2 

The higher-order integrals rv- 1+k ^ may therefore be written as: 



n 



(l+fe) _ 



1 + k 



(31) 
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The saturated effective number n s can then be reiated to the unsaturated effec- 
tive number n through an analytical expression: 

n s = nF (s m ) (32) 

F( Sm ) = l-^ + ^-^£ + .. = J-log(l + 2 Sm ) (33) 
2 3 4 2s m 

This may be compared with the corresponding equation in the plane-wave 
model: 

n s = nF pw (s) (34) 

F pw (s) = 1 - 2s + (2s) 2 - (2s) 3 + .. = — l — (35) 

The two expressions nearly coincide in the low saturation regime, where only 
the linear index and the Kerr term are appreciable, provided that an effective 
value s e of the saturation parameter is chosen in the plane-wave model such 
that: 



s m = 2s c (36) 

In more physical terms, this condition means that, in the Kerr limit, the plane- 
wave and Gaussian results have to be compared for equal mean intensities, i.e. 
equal intensities after averaging over the Gaussian profile. As soon as the satu- 
ration parameter is appreciable in contrast (s w 1 or s > 1), the two expressions 
differ. More strikingly, the expression ( ^3|) which accounts for transverse mode 
profile does not tend towards the plane-wave expression (|35| ) when the beam 
size parameter grows. In this regime also, valid conclusions about the effect of 
gaussian beam structure may be reached only after a criterium has been chosen 
to compare plane-wave and Gaussian-beam spectra. In the following, we will 
always compare spectra with equation (|3^) obeyed. 

The modification dA of the intracavity Gaussian mode amplitude due to one 
passage in the atomic medium may now be written: 

dA _ 3\ 2 nF(s m ) 

A 4ttS 1 + iS ( ' 

We easily obtain the steady state solution of the differential equation ( |l5| ) for 
the intracavity amplitude A: 



s /2~K>rcA m = KTc{l + i<p)A-dA (38) 
which may be written as: 



1 + 10 



(39) 
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Prn denotes the dimensionless parameter for the on-axis intracavity field (cf. 
equation (p^)) and f3 ln a dimensionless parameter for the input field defined in 
an analogous manner: 



/3 in = 



3A 2 A 1 ' 



(40) 



AttS 

C is the cooperativity parameter, which is proportional to the effective number 
n of atoms present in the probe beam: 



C 



Air ST 

T is the intensity transmission of the coupling mirror: 



(41) 



T = 2ktc 



(42) 



We then derive the relation between the input intensity and the intracavity 
intensity: 



^i/? in r=i/3mi 2 



(1 + i(j>) + 2C 



F (s m ) 



1 + iS 



(43) 



or: 



Y = X 



(1 + CF' {X)f + (cj> - SCF' (X)f 



(44) 



where we have replaced the field parameters by normalized intensities X and Y 
defined from the intensities averaged over the Gaussian profile: 



X= 1 -\f3 m \ 2 = l -s m (l + 5^) 



F> (X) 



2F(s m ) 1 / AX 

TTW = 2X log { 1 + TTP 



(45) 
(46) 
(47) 



9 E0 21, 231. It must be 



Equation (E3) is identical with results of [1 
compared with the corresponding equation in the plane-wave model: 



Y = X 



(1 + CF' (X)) 2 + (4>-5CF' (X)) 2 



(48) 



where the normalized intensities X and Y are now defined from the position- 
independent intensities in the plane-wave model: 



X= \[3 C \ 2 = s c (l + S 2 



(49) 
(50) 
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pi , x % _ 2f pw (fe 



pw 



< J > " " < 51 » 

As long as the saturation parameter remains small (X -C 1 + <5 2 ), the Insta- 
bility curves (Q) and (p|) are identical. Note that this regime corresponds to 
a low-saturation limit where the mean population of the atomic excited state 
defined in the plane- wave model remains small: 

A « 1 (52) 



1 + S 2 + 2X 

In the following, we will focus the attention on this regime where the predictions 
of the plane-wave and Gaussian-wave theories coincide for the mean fields. As 
already emphasized in the introduction, this does not imply that the predictions 
will also coincide for the fluctuations. 

Before coming to the study of fluctuations, we want to give conditions for 
obtaining a large squeezing. A first condition is that a bistability turning point 
is approached and references therein). This condition can be formulated in 
more quantitative terms by writing the slope of the bistability curve: 



rlV 

— = (l + a) 2 + (0-a<5) 2 -C 2 (53) 



where: 



a = 



2C (1 + S 2 ) 
(1 + S 2 + 2X) 2 



(54) 



C = -i^M+ZL (55) 
(l + 5 2 + 2X) 2 

The coefficient a describes the damping of intracavity field due to atomic absorp- 
tion, normalized to the damping due to the coupling mirror while the coefficient 
£ represents the normalized strength of the Kerr coupling due to atomic non- 
linearity. The presence of a turning point (dY/dX — 0) requires that the Kerr 
coupling is large enough so that third term in the right-hand side of equation 
(H|) is larger than the first one: 

(>l + a (56) 

and then that the detuning <f> of empty cavity is chosen such as to reach one of 
the two turning points: 



4>-a6 = ±y( 2 -(l + a) 2 (57) 

The second condition for getting a large squeezing is that bistability is mainly 
in the dispersive regime, which means that the effect of atomic absorption on 
the field remains small: 



10 



a « 1 (58) 

If this condition is not satisfied, atoms efficiently couple vacuum fluctuations 
coming from the empty modes responsible for spontaneous emission into the 
intracavity field fluctuations, leading to a degradation of the expected squeezing. 

In the low-saturation regime (X <C 1 + S 2 ), conditions ( |58"| , j56"| ) imply a large 
atom-laser detuning (<5 2 3> 1) and they may be rewritten: 

5 < 2C < S 2 (59) 

£ < X (60) 

In the following, we will consider that these conditions delineate the region 
where a large squeezing may be obtained. It is in particular worth noting that 
the intensity parameter X lies in the domain 5 <C X <C S 2 . 



4 Description of fluctuations in a Gaussian- wave 
theory 

We are now going to study the quantum fluctuations of the field. We will closely 
follow the reference fl6f where quantum fluctuations were derived from a linear 
input-output theory with a plane-wave model. As in the previous analysis of 
mean values, the reflection-transmission equations for the field fluctuations on 
the coupling mirror are the same as in the plane-wave model. In contrast, the 
effect of the atoms on these fluctuations is modified, due to the non-linearity of 
the interaction and to the Gaussian profile of laser intensity. 

In order to discuss the transformation of field fluctuations, we now introduce 
a new notation for the positive-frequency and negative-frequency components 
E (r,t) and E^ (r,t) of the field operator (see equation ([!])) which are coupled 
by the nonlinear interaction. They will be considered as the components £ a of 
a two- fold column matrix £: 



£ (r,i) = 



E (r,i) 
£t ( r ,t) 



(61) 



We will also denote £ a = ft the components of the adjoint line matrix £\ 

The field components are written in terms of mode amplitudes and Gaussian 
mode functions: 



£ a (r,t) = J^A a (t)U a (v) (62) 

The two components U a (r) coincide with the functions u (r) and u (r)* defined 
previously in equation (f 
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U a {v) 

The same notation will be used also for the positive-frequency and negative- 
frequency mode amplitudes: 

(65) 

as well as for the components of the atomic dipole measured as a dimensionless 
number: 

(66) 

(S (t) and (i) have the same definition as in |l6| |; they will be denoted S a (t) 
in the following). 

As for the mean values, the effect of an optically thin atomic layer of length 
dx may be described by a local modification of the field fluctuations: 

58 a (r, t) -> 88 a (r, t) + d (6£ a (r, *)) (67) 

We have introduced here a generic notation for the fluctuations SO of an oper- 
ator O: 




exp (-ie a tp (r)) 



(63) 
(64) 



A(t) 



A{t) 
At ( t ) 



S(t) 



Sit) 
St (t) 



SO = 0- (O) (68) 

The field modification d (S£ a ) may be separated into two parts representing 
respectively the linear response d (S8^j of the atomic dipoles to the intracavity 
field fluctuations and the spontaneous emission noise d(S8^), that is the linear 
response of atomic dipoles to the vacuum fluctuations in the field modes coupled 
to the atoms which are different from the intracavity modes (2^, |2^| : 

d (S8 a (r, <)) = d (88* (r, t)) + d (58 s a e (r, t)) (69) 

The spontaneous emission term is uncorrelated with the intracavity field fluctu- 
ations. Note that we disregard the effect of the higher-order intracavity modes, 
in consistency with our assumptions on the treatment of transverse mode struc- 
ture. We now evaluate these two terms. 

The first term may be written in terms of linear susceptibility functions Xa : 

d « r (r, <)) =—^rP (r) dx J Y^xi (r, r) S£p (r, t-r)dr (70) 
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The susceptibility function Xa ( r > T ) is deduced from commutators of the atomic 
dipoles evaluated at different times: 

xi(r,r)=ie(r)([S a (r),S^(0)]) r (71) 

9 (t) is the Heaviside function. The symbol (...) r means that the correlation 
functions are evaluated from the solution of the optical Bloch equations for an 
atom located at point r. As usually, the optical Bloch equations are solved in 
the frequency domain rather than in the time domain. We will give below the 
solution for the Fourier transforms Xa ( r > M) °f Xa ( r J r ) defined according to 
the general prescription: 

/(*)=/ ^e^/[ U ] 

The field modification d will then be written in the frequency domain: 

d « r (r, [u])) =~^-P (r) dxJ2x f i (r, [w]) 5E fi (r, [w]) (72) 

Clearly, the susceptibility function will depend upon the position of the atom 
because of the Gaussian mode structure of the field, which will make a difference 
with the plane- wave model of reference ]l6fl . 

Before considering this problem of position-dependence of the susceptibility 
function, we come to characterize the second fluctuating term d (SS^) appearing 
in equation (|69|), which represents the spontaneous emission noise. Assuming 
a dilute atomic medium such that the fluctuations of the different dipoles are 
independent of each other, we obtain field fluctuations by summing up the 
contributions of all dipoles: 



(d (5£? (r, H)) d (6£V sc (r', [u/]))) = 2n5 ( w + «') 6 (r - r') — p (r) dx (r, [w]) 

(73) 

The correlation functions erf describe the fluctuations of the dipole components 
evaluated for one atom at point r: 

a£(r,T) = (tfS a (r)W(0)> r (74) 

The analytical expressions of the Fourier transforms erf (r, [ui] ) , as well as those 
of Xa ( r > can be obtained in terms of a single atomic function which 
will be given in the next section. 

As already carried out for the mean values, we now solve the problem of 
the transverse dependence of correlation functions by supposing that the fluc- 
tuations radiated by the atoms into the higher-order modes have a negligible 
effect on intracavity field fluctuations, and furthermore are not matched to the 
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detection setup. We will therefore obtain the modification d (SA a ) of the ampli- 
tude fluctuations due to the atomic medium by projecting the field fluctuations 
d{5£ a ) onto the Gaussian mode functions: 

d(SA a (<)) = J d(5£ a (r,t))U*(r)dy dz 

Considering first the linear response term, we project the effect of a thin atomic 
layer of length dx (cf. equation (j70|)) and sum up over the whole atomic medium. 
We thus write d {8A 1 ^) in the frequency domain in terms of a new susceptibility 
function \a M : 

Q\2 

d {8Al [w])=- ^n^xi M SAp [u] (75) 



which is defined as the local susceptibility function averaged over the beam 
profile: 

Xl \u\ = -\ xi (r, [u])p (r) / (r) e^^dr (76) 



For the spontaneous emission noise, we once more disregard the fluctuations 
radiated by the atoms into the higher-order modes, then project the field fluc- 
tuations d5£^ onto the Gaussian mode, and finally sum up the independent 
contributions of all atoms. We thus obtain the contribution of spontaneous 
emission noise in terms of a correlation function er^ [lj] averaged over the trans- 
verse profile of the beam: 

(d (6A% C H) d (SA 13 sc [u/])) = 2nd (u + c/) — nai [u] (77) 

ai H = - f o% (r, [u])p (r) / (r) e «'«-' l >M) dr ( 78 ) 
n J 

Now, the computation of quantum fluctuations in the output field proceeds 
along the same lines as in the plane- wave model |]l6f , provided that the following 
modifications are accounted for. First, field fluctuations are understood as refer- 
ing to the Gaussian mode amplitudes, for the intracavity field 5A a as well as 
for input and output fields 8 A™ and <5„4° ut . Then, the plane- wave susceptibility 
function \a an d noise spectrum are replaced by the averaged expressions Xa 
and IfP given respectively by equations (76) and (fftj). This leads to a simple 
and natural generalisation of the results of the plane-wave model. When trans- 
verse beam structure is accounted for, the position-dependent expressions for 
atomic susceptibility function and atomic noise spectrum have to be averaged 
over the beam intensity profile. The same property was in fact used in pub- 
lished calculations accounting for Gaussian profile of the laser beam ||l], p3f : 
in these papers where quantum stochastic methods were used, the squeezing 
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spectra were obtained by averaging the drift and diffusion matrices which de- 
scribe in these methods the effect of atoms. In the present paper, the fact that 
the contributions of elementary volumes have to be added independently has 
not been postulated. We have shown that this averaging property follows from 
a simple treatment of Gaussian profile where fields radiated by the atoms are 
projected onto the fundamental Gaussian mode. 

5 Evaluation of the averaged atomic spectra 

Before coming to the estimation of the squeezing spectra, we now give the 
expressions of the atomic spectra and deduce the expressions averaged over the 
intensity profile. We then compare the atomic spectra thus obtained with those 
derived in the plane-wave model. 

The local susceptibility function and noise spectrum Xa ( r , M) and (r, [u]) 
are obtained by solving optical Bloch equations for an atom located at point r. 
They have the same form as in the plane- wave model [TsJ , with a position- 
dependent field parameter /3 (r) however, and may be written in terms of the 
phase (f (r) of the Gaussian mode and of a function which only depends on 
the local intensity parameter: 



X^(r,M) = ^{Gt (r,[c])-G^(r,[- W ])*)e-^-^M r ) (79) 
af(r,M) = i(Gf(r,H) + G;(r,H)>- l(e "^ )p(r) (80) 




(81) 



where: 




iu>) (2 — ilu) 



+2|/3| 4 (2-^) 2 +4|/3| 6 
i(|/?| 2 ,M) = 2\(3\ 2 {l + S 2 + 2\P\ 2 )(2- l uJ)(l + tS) 



(82) 



+ |/3| 2 (1 + iS) 2 (1 + iS- ioJ) (2 - iuJ) 
+2\/3\Uu(l + t6)+4\l3\ 6 



(83) 




• 2 (|/?| 2 ,M) = \f3\ 2 (l-iS) 2 (l-iS-iu)(2-iu) 
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+2\f3\ 4 iZJ(l-i5)+4\f3\ 6 (84) 

p|(|/f,M) = 2|/?| 4 (2-z^) 2 +4|/3| 6 (85) 

g(|/3| 2 ,M) = (l + + 2 |/3| 2 ) 2 (2 - iZZJ) ((1 - ^) 2 + 5 2 ) 

+ (l + <5 2 + 2|/?| 2 ) 2 4|/?| 2 (l^) (86) 



S and (3 are the dimensionless parameters already denned for measuring laser- 
atom detuning and Rabi frequency respectively; ZJ is the noise frequency also 
normalized to the decay rate of the dipole: 



In order to perform the average over the beam profile, we first note that the 
phases do not enter the averaged expressions (compare equations (^,^||) with 
equations ( |79| , p0| ) ) . We then proceed by analogy with the computation of the 
factors rS 1 ^' which appeared in the expression (^?j) of the effective number of 
atoms, by developing the functions in terms of the intensity parameter: 

oo 

Gf(r,[,])^GfWH|i3(r)| a (88) 

Writing the position-dependent local intensity as in the previous discussion of 
the bistability curve: 

|/3(r)| 2 = |/3 m | 2 /(r) = 2X/(r) (89) 

we obtain: 



where: 




The plane-wave expressions \a M an d a a M written in equations ( |79| , jS0| ) 
may be identified with equations ( |79| , |80| ) with replaced by unity (global phase 
factors are unessential and may be disregarded). The definitions of X, standing 
for the position-independent intensity parameter |/3| 2 in the plane-wave model 
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and for its spatial average over the beam profile in the Gaussian-wave theory, 
have been precisely chosen so that the coefficients £o and £1 are equal to unity. 
This implies that the differences between the plane-wave and Gaussian-beam 
predictions for squeezing only come from the fact that the higher-order coeffi- 
cients (with k = 2, 3 . . .) differ from unity. The discussion of these differences 
may therefore be restricted to the analysis of the contributions to atomic spectra 
proportional to X 2 , X 3 . . .: the plane- wave predictions cannot be considered as 
reliable as soon as these contributions become appreciable. 

It could be expected that the contributions proportional to X 2 , X 3 ... to 
the atomic spectra become appreciable, as for the mean values, only when the 
excited state becomes significantly populated, that is when the intensity param- 
eter X approaches S 2 . The plane- wave predictions would thus be reliable in the 
whole range 5 < X < S 2 where a large squeezing is expected (see the discussion 
at the end of section 3). We now show that this is not the case. 

We first consider the limiting case, that we will call the Kerr domain in the 
following, where only the terms proportional to X° and X 1 are appreciable. 
In this case, the coincidence between the atomic spectra obtained in the two 
models is guaranteed. The linear terms (proportional to X ) are given by: 

G{ {0) = : .] ._ (93) 
1 1 + id - S3 

G M°) = G? (0) =G^ (0) =0 (94) 
and the Kerr terms (proportional to X 1 ) by: 

= ^t^L (95) 



(1 + 5 2 ) (X + iS-iuy 
1 



(1 - %S) ((1 - S3) 2 + <5 2 ) (1 - iS f (1-iS- S3) 
= „ , ..a' . f (97) 



(1 + i5 f (1 + iS - iuJ) 



G\ (1) = (98) 

These terms have a simple physical interpretation in terms of a diagrammatic 
perturbative expansion of the field-matter interaction (3^, [HJ . The linear terms 
are described by two-photon elastic scattering amplitudes which are resonant at 
frequencies close to the atomic eigenfrequency (J3 ~ 6) while the Kerr terms are 
described by four-photon inelastic scattering amplitudes which are resonant at 
frequencies either close to the atomic eigenfrequency (~ ~ 5) or symmetric to 
the atomic eigenfrequency with respect to the laser frequency (to ~ — S). The 
interpretation of squeezing in terms of these photon scattering amplitudes was 
already performed for the case of atoms coupled to field in free space [E6l Blf. 
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The previous equations give a more complete description of the effect of these 
amplitudes which may be used also for interpreting squeezing by atoms coupled 
to field in a cavity. An important consequence of this discussion is that the 
squeezing processing is resonant on the sidebands of fluorescence triplet (ul ~ 
±<5) and not on the central component (ZU ~ 0). In contrast, squeezing via 
bistability is usually studied for frequencies close to the laser frequency. This 
is particularly clear in the good-cavity limit where the entire cavity bandwidth 
is contained in the central component of the emission triplet. In this case, the 
previous equations have to be evaluated at zero frequency (u7 = 0) where they 
correspond to non-resonant processes. 

We come now to the analysis of the higher-order contributions to atomic 
spectra proportional to X 2 , X 3 . . ., and focus the discussion on the same regime 
where squeezing is studied for frequencies close to the laser frequency. In this 
regime, it is easily checked on the expressions of the atomic susceptibility func- 
tions Xa that the higher-order contributions remain negligible as long as X <C 8 2 . 
This is related to the general property of linear response theory which states 
that static susceptibilities evaluated at ZD = are directly related to differenti- 
ated forms of the relations between mean fields, and to the already known fact 
that the plane-wave model provides reliable expressions for mean fields when 
X -C S 2 . However, the higher-order contributions to the noise spectra erf may 
become appreciable even when 1 < i 2 . This can be shown by considering the 
atomic spectra at zero frequency (obtained by setting = in the unapproxi- 
mated forms (|Bl))): 



(l + 5 2 + 2 |/3| 2 ) (l + 5 2 + 2\p\ 2 y 
Gl[Q] = ^21^(1 + ^^1^(1 + ^ + 2^ (100) 



l + <5 2 + 2|/?| 2 ) (l + £ 2 



•2 



GlM = Id' (!-*)' + * If (101) 



(l + S 2 +2\f3\- 
4|/?| 4 + 2|/3| 



Gi[0] = - '"' x3 (102) 

l + <5 2 + 2|/3| 2 

The evaluation of the susceptibility spectra relies only upon the first terms 
appearing in G\ [0] and G\ [0], so that the higher-order contributions to these 
terms are negligible when X <C S 2 . In contrast, the other terms, which enter 
in the evaluation of the noise spectra, may be appreciable even for X <C S 2 . 
Restricting our interest to this low-saturation limit, we may identify the main 
contribution to the squeezing process as the parametric terms proportional to 
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X/S 3 which appear in G\ [0]. Among the dominant noise contributions, there 
are the terms proportional to l/<5 2 in G\ [0] and X/S 4 in G\ [0] or G\ [0], but 
also the higher-order terms proportional to X 2 /S 6 or X 3 /5 6 which appear in the 
four functions G^ [0]. It is thus clear that these terms may not be disregarded 
for X ^ 5 2 . In particular, the noise term proportional to X 3 /S 6 reaches the 
magnitude of the parametric Kerr coefficient X/S 3 for X S 3 ^ 2 . 

This noise term is not accounted for in a Kerr model and is found to be 
responsible for excess noise and degradation of the expected squeezing [|| . From 
the qualitative evaluations of previous paragraph, one deduces that the noise 
correlation functions have to be suspected in the domain S < X < 5 2 , where 
a large squeezing is expected from plane-wave computations (see conditions 
( p9 60)). This does not imply that plane- wave predictions of squeezing are 
wrong in the whole domain, but this forces us to perform the Gaussian-wave 
computations in order to get reliable expectations. In the next section, we give 
a quantitative estimate of squeezing degradation due to this effect and delineate 
the domain of validity of plane- wave computations. 

The conclusion of this discussion is that, although the plane- wave model 
provides reliable estimations for the bistability curve and for the susceptibility 
functions in the whole low-saturation domain X <C 5 2 , this is not the case for 
the noise correlation functions. A qualitative interpretation of this result may 
be found by coming back to photon scattering amplitudes. As already discussed, 
the two-photon and four-photon amplitudes are resonant on the sidebands of the 
fluorescence triplet. In contrast, there exists higher-order scattering amplitudes 
which give rise to inelastic fluorescence on the central component of the triplet 
[ p0| . Although these amplitudes are proportional to higher-order powers of 
X, they are favored with respect to the lowest-order amplitudes, due to their 
resonant enhancement around the central frequency. This explains why they 
can have an influence in the low-saturation domain. 



6 Squeezing spectra 

We come now to the computation of quantum fluctuations in the output field. 
As already discussed, this computation proceeds along the same lines in the 
Gaussian-wave computations as in the plane-wave model, provided that the 
matrices Xa an d &a are replaced by the averaged expressions Xa an d &a ■ 

The covariance functions for the input or output fields are defined in terms 
of noise matrices: 

(5A% [u] 8A P in [a/]) = 2^5 (co + Vf in M (103) 

(5A° a ut H 6A P out [u/]) = 2ir5 (w + u') Vf out [u] (104) 

Vf in and Vf out are the elements of square matrices V ln and V out which are con- 
nected through noise processing relations (equations (47-51) of reference pq|): 
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v out M = (fi M - 1) .v in M . (fi [u] - i) f + fi H .v at M ./i M f (105) 

(/i [w] — 1) and fj, [co] are transfer matrices which characterize noise processing 
by the cavity containing atoms and which depends on the susceptibility matrix 

X- 

M M = —. - .J - „ 1 (106) 

V at describes the noise added by the atomic medium, and which depends on the 
atomic noise matrix a: 

V at [u] = jCe.a [uj] .e (107) 
In these expressions, e is a diagonal matrix containing the coefficients e a : 



£ = 



1 
-1 



(108) 



It is worth stressing at this stage that, while the atomic spectra x an d <? are 
obtained by averaging the plane-wave spectra over the intensity profile, this is 
not the case for the noise spectra of the output fields. 

We may now give a quantitative estimation of the influence of the Gaus- 
sian profile upon the noise spectra for the output fields. The transfer matrices 
are not modified in the low-saturation limit, since they depend only upon the 
susceptibility matrix x- It follows that the first term in the right-hand side of 
equation ( |105| ) is not influenced by the Gaussian profile. In contrast, the second 
term, which is proportional to the atomic noise matrix a, has a larger degrading 
effect on squeezing in the Gaussian-wave theory than in the plane-wave model. 
Precisely, the squeezing degradation due to higher-order noise terms is greater 
in Gaussian-wave theory than in the plane-wave model; for terms proportional 
to X 2 (respectively X 3 ), this degradation is multiplied by £2 = 4/3 (respec- 
tively by £3 = 2). In the Kerr domain, those higher-order noise terms have a 
negligible influence on squeezing, so that the predictions of plane-wave model 
coincide with those of Gaussian-wave theory. 

A more detailed discussion may be given by studying the example of the 
optimum squeezing expected at zero frequency (oj = 0) on a bistability turning 
point (defined by equations ( [56|j5l| )). The optimum squeezing (variance of the 
quadrature component with minimal fluctuations, evaluated at zero frequency 
on a bistability turning point) reaches in the Kerr domain the value: 

v 1 + a 

This value is small, of the order of 2C/S 2 , when condition ( |59| ) is obeyed, that 
is when the effect of atomic absorption on the field remains small. Note that 
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a small value of S Q pt means an efficient squeezing of input field fluctuations by 
the bistable cavity. Assuming that the intensity parameter X lies in the domain 
S -C X <C S 2 , a close inspection of the equations describing noise processing 
by the bistable cavity shows that the optimum squeezing S op t accounting for 
higher-order noise terms in the plane-wave model may be written: 



oKerr i o/-^ 
'-' opt + z J ^ 

1 + 2JV 



Sopt — — ~ ; n j. x 2& (HO) 



The numerical factor / has the following value, of the order of unity (<j) is defined 
by®): 

f (111) 
1 + a 

The degradation due to higher-order noise terms is thus found to become ap- 
preciable when 2X 2 /S 3 ~ 2C/8 2 , that is when: 

X~VCS (112) 



The degradation of squeezing appearing in equation ( |110| ) is enlarged when 
passing from the plane-wave to the Gaussian-wave theory. It turns out that, 
in the limiting case where degradation is small, it may be attributed mainly to 
higher-order atomic noise terms proportional to X 3 , so that squeezing degra- 
dation is roughly multiplied by £3 = 2 due to Gaussian transverse profile. In 
other words, the plane-wave model systematically underestimates the degrading 
effect of higher-order contributions to atomic noise and does not provide reliable 
estimations for the squeezing when such contributions cannot be neglected. 

In contrast, we may consider now the condition which delineates the Kerr 
region: 

X^VCS (113) 

When this condition is satisfied, the optimum squeezing is not degraded by 
higher-order noise terms, so that the plane- wave computation provides a reliable 
estimation for squeezing. It is worth recalling here that this condition has been 
derived in the domain S <C X <C S 2 where a large squeezing is expected. In the 
domain of very low intensity X < 5, different conclusions would be drawn for 
reliability of the plane-wave model, but a poor squeezing would be predicted. 

The previous discussion was restricted to the particular case of zero fre- 
quency where analytical expressions may be handled more easily. Figure 1 
shows spectra for the optimum squeezing evaluated as a function of noise fre- 
quency u>. The three spectra correspond to the same value for the atomic ab- 
sorption coefficient a and to the same working point on the bistability curve; we 
have chosen the bistability threshold such that there is only one turning point, 
which corresponds to C — 1 + a and (ft = a8 (/ = 1 in this particular case). 
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Figure 1: Spectra for the optimum squeezing as a function of normalized noise 
frequency uj = w/7, with a decay constant for the intracavity field k = IO7; the 
three spectra correspond to the same atomic absorption coefficient (a = 0.25) 
and to the same turning point at bistability threshold (£ = 1 + a; cf> = ad); 
the spectrum (i) is computed for a large detuning lying in the Kerr domain 
(S = 1000; the other parameters deduced from the values of 6, a and ( are 
X = 2500, C = 126253, 4> = 250); the two other spectra are computed for a 
lower detuning lying outside the Kerr domain (6 = 100; X = 250, C = 1378, 
4> = 25) and result from plane-wave computations (ii) and Gaussian-wave com- 
putations (iii). 
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The first spectrum (trace (i)) is computed for a detuning sufficiently large so 
that the Kerr limit is reached. In this case, the spectra computed from plane- 
wave model and Gaussian-wave theory cannot be distinguished. The two other 
spectra are computed for a lower detuning, such that higher-order noise terms 
degrade optimum squeezing. The results of plane- wave computations (trace (ii)) 
and Gaussian-wave computations (trace (hi)) now differ, revealing in particular 
a larger degradation when transverse structure is accounted for. In addition, 
these spectra clearly show, in consistency with the discussion at the end of sec- 
tion 5, that the excess noise responsible for squeezing degradation is associated 
with the central peak of the atomic fluorescence spectrum, whose width equals 
27, where 7 is the decay constant of the atomic dipole, in the limiting case 
of low saturation |32| . Note that the spectra of Figure 1 have been computed 
in the bad-cavity limit where the decay constant k of the intracavity field is 
larger than 7. In the opposite good-cavity limit, the resonant behaviour of the 
higher-order excess noise would not be apparent on the squeezing spectrum. 



7 Conclusion 

We have used the linear input-output formalism generalized to incorporate the 
treatment of atomic fluctuations to discuss the effect of Gausian transverse 
structure on squeezing with two-level atoms. We have shown that this effect is 
described quite naturally as a spatial averaging of the susceptibility and noise 
functions which appear in the linear input-output formalism. This property 
has not been postulated a priori, but derived as a consequence of the simple 
assumptions made in the treatment of transverse profile: higher-order transverse 
modes have been supposed either far from resonance, or of such a high order 
that they are not efficiently coupled to the Gaussian mode. 

The mean fields and the bistability curves deduced from plane-wave and 
Gaussian-wave calculations coincide in the low-saturation limit where the pop- 
ulation of the atomic excited state remains small (X <§; S 2 ). We have shown 
that this is also the case for the atomic susceptibility functions evaluated around 
zero frequency. This is related to the general property of linear response theory 
that the static susceptibilities are differentiated forms of the relations between 
mean fields. In contrast, the limit of low saturation is not sufficient for ensuring 
the coincidence of plane-wave and Gaussian-wave predictions for atomic noise 
functions. This is due to the existence of higher-order noise terms, arising from 
the central peak of the fluorescence spectrum. These terms are able to degrade 
the optimum squeezing when X reaches the value V C5. We have shown that 
these terms have a larger degrading effect in the Gaussian-wave calculations 
as in the plane- wave ones. This implies that plane- wave predictions of a large 
squeezing are reliable only in the Kerr domain, that is for X <C V CS. This 
criterium for reliability of plane- wave model appears quite different from those 
published in previous references (2ll E2l E3j . 
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